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Abstract 

We investigate the feasible limits for realising a continuously evaporated atom laser with high- 
temperature sources. A plausible scheme for realising a truly continuous atom laser is to outcouple 
atoms from a partially condensed Bose gas, whilst continuously reloading the system with non- 
condensed thermal atoms and performing evaporative cooling. Here we use quantum kinetic theory 
to model this system and estimate feasible limits for the operation of such a scheme. For sufficiently 
high temperatures, the figure of merit for the source is shown to be the phase-space flux. The 
dominant process limiting the usage of sources with low phase-space flux is the three-body loss 
of the condensed gas. We conclude that certain double-magneto-optical trap (MOT) sources may 
produce substantial mean condensate numbers through continuous evaporation, and provide an 
atom laser source with a narrow linewidth and reasonable flux. 
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I. INTRODUCTION 



The experimental realisation of Bose-Einstein condensation (BEC) in dilute alkali gases 
PP has opened up the study of the atom laser: a coherent beam of atoms analogous to the 
optical laser. Like optical lasers, it is hoped that these atomic sources will demonstrate mode 
selectivity and high spectral flux [2j. The simplest method of producing a spatially coherent 
atomic beam is to couple atoms out of a trapped BEC [3]. Rapid out-coupling of a BEC 
forms a coherent atomic beam with a spread in momentum as large as the trapped BEC, 
whereas coupling the atoms out more slowly reduces the output linewidth at the expense of 
reducing the overall flux j3]. These atom lasers are equivalent to Q-switched optical lasers, 
which do not exhibit gain-narrowing. Continuously pumping the atom laser to produce a 
stable output provides the obvious benefit of higher flux, but may also improve stability and 
linewidth of the output beam. In a gain-narrowed optical laser, a higher pumping rate both 
increases the total flux and reduces the linewidth of the output, producing a dramatically 
increased spectral flux [5]. An atom laser with gain-narrowing must have a saturable, Bose- 
enhanced pumping mechanism that operates simultaneously with the damping [2J. This 
paper examines the limits on thermal sources that can be used to produce a continuous 
atom laser through the process of continuous evaporation. 

The two essential steps towards providing a continuous pumping mechanism for an atom 
laser are: (1) the provision of atoms from an external source to the atomic trap, and (2) a 
process that causes at least some of those atoms to make an irreversible, stimulated tran- 
sition into the BEC. Continuous delivery of ultracold atoms has been demonstrated in a 
number of experiments [HH9], and is an important component of thermal atomic interferom- 
etry. Sequential reloading of a target BEC was achieved using optical tweezers j6], where 
a series of source condensates were joined by manipulating the trapping potentials, and 
excitations were subsequently removed by further evaporative cooling. This milestone ex- 
periment maintained the condensate fraction, and therefore the potential flux of a potential 
atom laser. However, an atom laser produced from such an experiment would not possess 
the desired narrow linewidth, as merging two coherent sources by manipulating the poten- 
tial is not a Bose-enhanced process. This means that the phase of the lasing mode would 
diffuse at a rate proportional to the replenishment rate. In this paper it is shown that under 
certain conditions, a similar experiment using an ultra-cold thermal source would be able to 
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pump the target BEC and maintain a significant BEC population using a phase-preserving 
Bose-enhanced scattering process. 

Atom-stimulated transitions into the condensate can be made irreversible by coupling 
to a reservoir. There are two possible reservoirs: the empty modes of the electromagnetic 
field accessible via a transition from an excited atomic state (as used in optical cooling), 
or the empty modes of the atomic field (accessible via evaporation). Optical cooling is an 
obvious candidate for a continuous process, but it has not yet achieved BEC. This is because 
condensates are very sensitive to the resonant light that is emitted by the spontaneous 
emission step in the cooling, as the recoil energy of a single photon is significantly larger 
than the energy per particle. Several early atom laser models were based on optical cooling 
from a thermal source pU] , but they have not been experimentally realised. The only optical 
cooling method that has led to an increase in BEC number required a pre-condensed source 
in the same trap 

The last stage of cooling for production of all current BECs from thermal sources has 
been mediated by atom-atom interactions. This has the advantage that it can be performed 
without the presence of resonant light, but the obvious disadvantage that it relies on the 
system approaching thermal equilibrium, and will therefore be reversed by the addition of 
atoms above the condensation temperature. However, if the atoms can be supplied at the 
same time as the system undergoes a forced evaporation process, then even if the source 
atoms are above the condensation temperature, it is possible under some conditions to 
produce a net gain for the condensate mode. It is therefore possible to make evaporative 
cooling operate in a continuous fashion. The focus of this investigation is to determine 
exactly which existing thermal sources might be feasibly used for such an experiment. 

The quantum kinetic theory of dilute gas BEC has been successfully used to model 
condensate growth [121 [13]. We extend this model to include loss due to three-body inelastic 
scattering, which turns out to be a key determinant of the required phase-space flux required 
in the thermal source. In Sec. [TTJ we shall describe the schematic of our proposed pumping 
scheme, and the details of the quantum kinetic theory model. The results of the model are 
examined in Sec. IIHI 
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FIG. 1: Schematic of the experimental setup. 
II. SCHEME AND KINETIC THEORY MODEL 

The proposed scheme for a pumped atom laser is illustrated in Fig. [TJ It is very similar 
to the methods used to evaporate a thermal cloud to condensation in a magnetic trap and 
produce a (quasi- continuous) atom laser. The gain process for the condensate is the same 
Bose-enhanced scattering between thermal atoms and the condensate that drives condensate 
growth when evaporating to produce BEC [121 HH US]- This process becomes irreversible 
when one of the scattered atoms has enough energy to cross the evaporation surface and 
be removed from the thermal cloud. The loss of atoms from the thermal cloud is balanced 
by a replenishment process that couples the thermal cloud to a source of atoms at finite 
temperature. 

The atom laser beam itself is produced by large momentum-transfer Raman outcoupling 
from the condensate. Raman outcoupling improves the spatial properties of the beam [16J, 
and allows minimal outcoupling from the thermal cloud when the two lasers are focussed to 
intersect only in the immediate vicinity of the condensate. 

A non-equilibrium steady state will be reached when the rate of atom loss from the 
condensate due to outcoupling balances the rate of atoms gained due to growth collisions 
within the thermal cloud. If the evaporative surface is tuned so that atoms of energy e cnt 
and higher are rapidly and continually removed from the trap, then all collisions that result 
in an atom with energy greater than e cut will become irreversible. As e cut is lowered, a larger 
fraction of the scattering processes that leave atoms in the condensate mode will become 
irreversible. This suggests that there must be some value of e cnt for which the condensate 
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experiences net gain. What is not clear is whether the net gain can proceed efficiently, i.e. on 
a timescale much shorter than other losses from the condensate. Lowering e cut also reduces 
the total number of thermal atoms present in the steady state. In the limit that e cut reaches 
the condensate energy, there will be no background gas at all, and the condensate cannot 
experience net gain. We therefore expect that for a given set of parameters, there will be 
an optimal value for e cnt that maximises the net gain. This value can be determined by a 
model based on quantum kinetic theory (QKT) [T3] . 

Our model is an extension of the kinetic models described in [121 E] , which was success- 
fully used to study condensate growth in an experiment in which a cloud of thermal atoms 
just above condensation temperature were shock-cooled through the BEC transition j22j |23] . 
After cooling, the atoms were left to equilibrate, with condensate formation being driven by 
the same collisional processes that would drive condensate growth in the proposed pumped 
atom laser experiment described in the previous section. To fully describe this proposed 
experiment, the kinetic model must be modified to include the effects of the replenishment 
and outcoupling processes, as well as the loss due to three-body recombination. 

Three-body recombination is the dominant loss process in typical BEC experiments [18] , 
and it provides the main constraint on the size of the lasing mode. In fact, if the outcou- 
pling from the condensate were the only non-negligible loss process from the trap, the largest 
condensate would be formed in the absence of evaporation. Furthermore, the condensate 
number would be independent of the temperature of the replenishment source, depending 
only on the flux of atoms delivered to the system and the outcoupling rate from the con- 
densate. This unphysical result is due to fact that in the absence of a density-dependent 
loss process, simply increasing the density is a feasible method of approaching degeneracy. 
Density-dependent loss puts a limit on the physical density, and thus results in the expected 
trade-off between temperature and flux density into the model. 

The starting point of the kinetic model is a separate treatment of the thermal and con- 
densed components of the system. The condensed component is assumed to be a quantum 
fluid obeying a Gross-Pitaevskii-type equation, however we make a further approximation 
and assume that the condensate is sufficiently occupied that it has a Thomas-Fermi pro- 
file [T9| Chapter 6]. The condensate dynamics are then fully described by the number of 
condensed atoms N (t). 

The thermal cloud is described by the Hartree-Fock approximation [T51 Chapter 8] , which 
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assumes it is comprised of particle-like excitations moving in the effective potential of the 
harmonic trap plus condensate mean field. To reduce the dimensionality of the full phase- 
space distribution function for the thermal cloud f(r,p,t), it is assumed that the system is 
ergodic, i.e. that all points in the phase space having the same energy have equal population 
[20] . Under this approximation, the thermal cloud is then described by its energy distribution 
function g(s, t) and the density of states p(e, t). The assumption of ergodicity has been shown 
in the past to give good agreement with experiment when asymmetric spatial or momentum 
dynamics are not significant [121 E] . Note that the time- dependence of the density of states 
p(e, t) comes from the contribution of the condensate mean field to the effective potential 
experienced by the thermal atoms. 

As the model presented here is very similar to that presented in [12] with some additional 
terms, a derivation of the common terms is omitted. As a summary, the derivation proceeds 
by taking a semiclassical Boltzmann equation for the phase-space distribution function of 
the thermal cloud f(r,p,t) including collisional terms and using the ergodic approximation 
to obtain an equation of motion for the energy distribution function g(e, t). This equation is 
self-consistently matched with a Gross-Pitaevskii equation for the condensate before making 
the Thomas-Fermi approximation to obtain an equation of motion for the number of con- 
densed atoms N (t). An example application of this method to derive the appropriate terms 
for three-body loss is given in section [3] of the appendix. Further, a detailed discussion of 
this theory is given in the review article [2"Tj . 

Separating the contributions of the different processes involved, the equations of motion 
for the model for a collision-driven pumped atom laser considered here are 
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where the subscripts 'thermal-thermal' and 'thermal-condensate' denote Bose-enhanced col- 
lisional processes between atoms in the corresponding states [Figs.|2](a) and (b), respectively], 



the subscript '3-body loss' indicates the contribution due to three-body recombination, the 
subscript 'replenishment' indicates the contribution due to the replenishment of the thermal 
cloud [Fig. |2[d)], the subscript 'outcoupling' indicates the contribution due to outcoupling 
from the condensate to form the atom laser [Fig. [2[e)], and the subscript 'redistribution' 
indicates the contribution due to the redistribution of population in energy space due to the 
changes of the energies of the occupied levels as the mean-field of the condensate changes 
[Fig. |2](c)]. It is assumed that atoms with energy greater than the evaporative energy cut-off 
e cut are removed from the system sufficiently quickly such that g(e > e cut ) = 0. 

The forms of the 'thermal-thermal', 'thermal-condensate' and 'redistribution' terms in 
([TJ are given in the appendix, and derivations are given in [12]. 

The outcoupling process from the condensate is modelled as a simple linear loss process 

with corresponding rate constant 7, 

dN 
~dT 

Modelling the outcoupling in this way neglects any outcoupling from thermal modes. This 
is a reasonable approximation if focused Raman lasers are used for the outcoupling, which 
only intersect in the immediate vicinity of the condensate. It also requires that the dominant 
losses from the thermal cloud are evaporation, coupling to the condensate and three-body 
loss. 

The thermal cloud is modelled as being continuously replenished from a source that 
provides a constant flux $ of atoms at a temperature T. To avoid tying the model to any 
particular replenishment mechanism, we assume a best-case scenario in which each energy 
level e in the source is coupled directly to the level in the thermal cloud with the same 
energy above the condensate chemical potential p(t), i.e. the lowest energy level of the 
source (e = 0) is coupled directly to the lowest energy level in the trap (e = pit)). This 
simple model gives the form of the contribution due to replenishment as 

d(p(e,t)g(e,t)) 



= ~lN . (2) 

outcoupling 



dt 



= Tp (e-p(t))g T (s-p(t)), (3) 

replenishment 

where po(e) is the density of states in the absence of a condensate, gr^) is the Bose-Einstein 
distribution at temperature T, and T is a rate constant such that 

/■oo 

r/ p (e)g T (e)de = $, (4) 
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FIG. 2: Schematic of processes involved in the evolution of the kinetic model described by ([I]). 
The upper shaded rectangle in each subfigure represents the energy distribution function g(e, t) of 
the thermal cloud, and the bottom dark blue rectangle represents the condensate with occupancy 
No(t) and chemical potential e = n{t). Figures (a) and (b) represent collisional processes involving 
two thermal atoms and one thermal and one condensate atom respectively. Figure (c) represents 
the change in the energy distribution function g(e, t) if the condensate occupation (and hence 
chemical potential) increases, raising the energies of every energy level. Figure (d) represents the 
replenishment of the thermal cloud from an atomic reservoir, and Figure (e) represents outcoupling 
from the condensate mode to produce the atom laser. 

where <3> is the flux of atoms from the source before evaporation. The derivation of the 
contributions to ([lj due to three body loss are given in the appendix. 

We summarise here the approximations made in obtaining the kinetic model Q: 

1. The energy scale of the thermal cloud is large enough that all excitations are particle- 
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like and not collective excitations such as phonons. Phonon-like excitations are only 
important for particle energies e < 2/i(£) [TjJl §8.3.1]. Hence, we require the that the 
energy scale for the thermal cloud e cut be much larger than 

2. The phase-space distribution of the thermal cloud is ergodic and hence is purely a 
function of energy. This assumption is true at equilibrium, however it needs some jus- 
tification when used in non-equilibrium scenarios. In this case, asymmetric behaviour 
of the condensate is not expected in either position or momentum space. 

3. The condensate density is sufficiently large that it is well-described by a Thomas- Fermi 
profile. This approximation is justified as it is only the large-condensate limit that is of 
interest as a large condensate will be necessary for the production of a high-flux atom 
laser in this scheme. In making this approximation, the effects of both the normal and 
anomalous densities of the thermal cloud on the condensate have also been neglected. 

4. Evaporation occurs on a time-scale faster than collisions. This is the usual require- 
ment during evaporation to condensation, and so should be satisfied in the proposed 
experiment. 

III. RESULTS 

For a given trap geometry, the model is fully defined by the flux of replenishment atoms $, 
the temperature T of those atoms, the energy of the evaporative cut e cut , and the outcoupling 
rate from the condensate 7. In this section the results of the kinetic model for some 'typical' 
parameter values are presented, and the dependence of the model on each of the parameters 
is examined. 

Our numerical simulations are based on a trap and conditions similar to that of [22], 
who precooled a cloud of 87 Rb atoms to an initial temperature slightly greater than the 
critical temperature before performing evaporative cooling to study condensate growth. The 
trap in the experiment was axially-symmetric with radial and axial trapping frequencies of 
bj r = 2n x 110 Hz and u z = 2tt x 14 Hz respectively. 

To solve the kinetic model numerically, Eq. ([T]) is discretised along the energy dimension 
and the resulting coupled differential equations are solved with an adaptive fourth-fifth 
Runge-Kutta [21] method. Our results are mainly concerned with the steady-state of the 
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FIG. 3: Results of the kinetic model for $ = 8.4 x 10 5 atoms/s, T = 540 nK, e cut = 3k B T & 610fiuJ, 
and 7 = 0.3 s . Figure (a) highlights the dynamics of the occupation of the thermal energy levels 
for t < 0.5 s, while (b) illustrates the equilibration of the total and condensed atom numbers over 
~ 10 s. The energy distribution at t = s is a truncated Bose-Einstein distribution containing 
(before truncation) N = 4.2 x 10 6 atoms at T = 540 nK. 

kinetic model, which we define as being reached when the condensate number has changed 
by less than 0.1% or 1 atom in 100 ms. The initial state for the simulation is chosen to 
be a truncated Bose-Einstein distribution containing (before truncation) iVi n itiai = 4.2 x 10 6 
atoms at the same temperature as the replenishment reservoir. This state is chosen as a 
representation of the steady-state of the system prior to evaporation. In the trap considered, 
the critical temperature for 4.2 x 10 6 atoms is T c = 400 nK. 



A. Typical results and parameter studies 

As a depiction of the 'typical' time-dependence of the results obtained from the kinetic 
theory model (111) , we consider the case of pumping the system continuously with a source 
such that the initial number N = 4.2 x 10 6 is transferred to the system once every 5 seconds 
giving a flux of $ = 8.4 x 10 5 atoms/s. The temperature of the replenishment source is 
chosen to be T = 540 nK, 60% above the condensation temperature of the system before 
evaporation. For the remaining model parameters, we choose the evaporative cut-off to be 
£ cu t = 3k B T, and the outcoupling rate from the condensate to be 7 = 0.3 s" 1 . 

Figure [3] illustrates the results of the simulation of this system. Fig.[3ta) shows the energy 
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distribution of the thermal cloud cooling from the initial truncated Bose-Einstein distribution 
to a distribution with a lower average energy per particle. Figure ^b) demonstrates that, 
despite pumping the system with an atomic reservoir above critical temperature, it is possible 
to reach a steady-state in which the condensate is macroscopically occupied. In this example, 
the steady-state condensate fraction is 33%. 

The non-equilibrium dynamics of the system are not the subject of investigation here, 
instead our interest is in the steady state itself, and in determining the feasibility of creating 
a pumped atom laser driven by a non-condensed atomic source. As a first step towards 
this investigation we consider the dependence of the equilibrium condensate number on 
the parameters of the system: T, e cut , and 7. The results of this parameter study are 
presented in Fig. [4j where we use the initial conditions of Fig. [3j 

In general the parameter dependences depicted in Fig. [4] are straightforward; adjusting 
each parameter causes a monotonic change in the equilibrium condensate number. Increasing 
the flux of atoms to the system increases the equilibrium condensate number [Fig. |4fa)], 
while increasing the temperature of the replenishment source or increasing the outcoupling 
rate reduces the equilibrium condensate number [Fig. |4](b) and (c) respectively]. The only 
behaviour that is not straightforward is displayed by Fig. |4^d) in which the dependence 
on the evaporative cut-off e cut is illustrated. For large e cnt , few atoms will be lost due to 
evaporation and the system will reach equilibrium when the flux of atoms into the system 
is balanced by three-body losses and outcoupling from the condensate. As e cut is reduced, 
more atoms are lost due to evaporation and the mean energy per particle reduces causing 
the condensate size to increase. As e cut continues to decrease, an increasing fraction of 
the replenishment atoms have an energy greater than e cut causing a lower effective atomic 
flux to be delivered to the system, hence reducing the steady state size of any condensate 
formed. These two competing effects are the origin of the existence of an optimum steady 
state condensate number as a function of e cut in Fig. |4^d). 

As discussed earlier, in the absence of three-body loss the equilibrium condensate number 
would continue to increase as e cut is increased, which would lead to the unphysical conclusion 
that evaporating reduces the equilibrium condensate number. This is demonstrated by the 
dashed line in Fig. |4^d) which asymptotes towards N = $/7 = 2.8 x 10 6 atoms in the 
limit e cut — > 00. As observed in the remaining panels of Fig. [4] (in which the effects of 
three-body loss have been included) three-body loss does not give rise to optimum values 
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FIG. 4: The dependence of the steady state condensate number iVo on the parameters of the 
quantum kinetic model Q. The steady state condensate number has a monotonic dependence on 
the replenishment flux <E> (a), the temperature T of the replenishment source (b) and the outcoupling 
rate 7 (c). For a given choice of the remaining parameters of the model there is an optimum e cn t 
(d) for which the equilibrium condensate number is a maximum. For each parameter being varied, 
the remaining parameters are chosen to be the same as for the results depicted in Fig. [3j The 
triangle in each plot marks the point that corresponds to the precise conditions of Fig. [3] in steady 
state. 

for the corresponding parameters of the model as it is only changes to e cut that affect the 
evaporative and three-body losses in contrary fashions. An increase in the replenishment flux 
will increase both evaporative and three-body losses. Similarly, changes to the temperature 
of the replenishment source or the outcoupling rate either increase or decrease both of the 
evaporative and three-body losses. 
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At this point it is clear that to have the largest equilibrium condensate number one should 
use a replenishment source with the highest possible flux and the lowest possible temper- 
ature. However any such recommendation would be disingenuous. It is an experimental 
reality that these parameters are not entirely uncoupled; while a 300 K oven might produce 
a significantly larger flux than a 50 mK 2D-MOT, it is simply not realistic to create a 50 
niK atomic source with the same flux as the 300 K oven. We will now examine this trade-off 
between the temperature and flux in the context of experimentally-realisable sources. 

B. Behaviour in the high-temperature limit 

In the previous section we investigated the dependence of the equilibrium condensate 
number on the model parameters. The physical question that we wish to address with this 
model is: What are the limits of a thermal atom source such that is can realise a pumped, 
continuous atom laser? 

Although it would be possible to create a pumped atom laser by combining condensates 
in a manner similar to the experiment by Chikkatur et al. |6J, such an atom laser would 
have significantly reduced phase-stability unless the replenishment process were essentially 
continuous. However, to replenish a condensate by collisional interactions with a continuous 
source of condensed atoms, the replenishment source would itself need to have many of the 
desired properties of a pumped atom laser. Instead, it would be preferable to use a source 
above condensation temperature for replenishment. 

We consider now the experimentally-relevant limit of replenishing the thermal cloud using 
a high-flux source of thermal atoms. For such sources two simplifications are possible. First, 
for temperatures greater than T c the Bose-Einstein energy distribution of the source gx{£) 
is well approximated by the Boltzmann distribution gr{£) ~ C e_/3e f° r some constant (, 
and (3 = (/cgT) -1 . Secondly, for high temperature sources the optimum evaporation cut- 
off e cut will be much smaller than the characteristic energy of the source fceT, and hence 
£ cu t <C ksT. From these simplifications it can be seen that the energy distribution below 
the evaporation cut-off is well described by the single parameter ( as gx{e < £ C ut) ~ 

At this point we have simply rewritten the temperature dependence of the replenishment 
source in terms of the parameter (. However, as the energy distribution of the replenishment 
source only affects the kinetic model through Eq. d3]), its influence on the system dynamics 
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is only through the combined quantity k = T(. An approximate expression for k in the 
Boltzmann limit directly in terms of relevant experimental quantities can be obtained using 
the definition Eq. Q, 



$ = T/ Po (e)g T (s)de (5) 

e 2 
2{tiuJ) 



r / HT^Ce-^de (6) 



oo 



r C^=^ / ^~^de (7) 



2{Wjf 

3 

rc (8) 



where u; = (u^u^Wz) 3 is the geometric mean of the trapping frequencies, and po(e) = 
is the density of states in a harmonic trap in the absence of a condensate [19]. 

We identify k as the phase-space flux of the source as it is directly related to the rate at 
which the phase-space density of the thermal source is delivered. For a trap of N thermal 
atoms at temperature T, the peak phase-space density w is [T9J Chapter 2] 

(£)'■ < 10 > 

If these iV atoms are delivered over a time r providing a flux $ = N/ r the peak phase-space 
flux is 

$ hr^ ( n ) 



The parameter £ is therefore the phase-space density ro. 

The phase-space flux k is a figure-of-merit for the thermal source. It quantifies the 
qualitative behaviour already known: for the same atomic flux $, a source with a lower 
temperature will result in a larger condensate [Fig. |4](b)]; and for the same temperature, 
a source with a higher atomic flux will also result in a larger condensate [Fig. |4j^a)]. The 
phase-space flux also describes exactly how a trade-off between the flux and temperature of 
the replenishment source will affect the equilibrium condensate number. If two sources with 
different fluxes and temperatures have the same value phase-space flux, then the equilibrium 
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condensate number produced by the two sources will be the same (assuming the high- 
temperature limit applies to both sources). Our interest is in determining what values of k 
are necessary to produce a pumped atom laser, and whether such values are achievable. 

For the limit of high-temperature atomic sources, we have reduced the four variables 
($, T, e cut , 7) required to define the model ([TJ down to three («, £ cut , 7). Of these three, 
our main interest is in the dependence of the system on the properties of the atomic source 
through k. In contrast, the dependence of the equilibrium condensate number on the out- 
coupling rate 7 is simple (see Fig. |4[c)) and the results would not be expected to change 
qualitatively with 7. It is therefore appropriate to choose a representative value for the 
outcoupling rate (here 7 = 0.3 s" 1 ) and focus on the remaining two quantities. 

As discussed in the previous section, there is an optimal choice for the evaporative cut-off 
Scut. Our interest here is in the best-case scenario: for a given thermal source, what is the 
largest condensate we can produce? To examine this question and to verify that k does fully 
describe the properties of the thermal source in the appropriate limit we have performed a 
parameter scan of the model Q for a range of fluxes 1.3 x 10 5 s _1 < $ < 5 x 10 10 s" 1 and 
temperatures 200 nK < T < 600/iK of the atomic source, for each combination determining 
the optimum evaporative cut e cut to give the largest steady-state condensate number. The 
results of this parameter scan are displayed in Fig. [5j 

The results illustrated in Fig. [5] are separated into three groups based on the ratio 
e cut / (k B T). The first group marked by black circles have e cnt /{k B T) < 0.1, and are the 
results for which the high-temperature limit can be considered to be a good approximation 
(satisfying the requirement e cut <C k B T). Here k completely determines the properties of 
the replenishment source. For this group of results, any equilibrium property of the system 
should appear to be a single (not necessarily straight) line when plotted as a function of k. 
The results in Fig. [5] demonstrate that these results can be viewed as a single function of k. 
The second group marked by red triangles have 0.1 < e CVLt /{k B T) < 0.5, and can be consid- 
ered to be the results for which the high-temperature limit is almost a good approximation. 
These results are reasonably close to the results of the first group, however there is a greater 
deviation for a given value of k indicating that the results can be almost seen as purely a 
function of k. All remaining results fall into the third group for which e cut /(k B T) > 0.5. It 
can be seen that these points correspond to a broad range of steady states, indicating that 
the replenishment source cannot be described by k alone. 
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FIG. 5: Steady-state condensate properties as a function of the phase-space flux k for the replen- 
ishment source. Figures (a) and (b) plot the steady-state condensate number on linear-linear and 
log-log scales. Figure (c) plots the equilibrium condensate fraction Nq/N. The results of the pa- 
rameter scan are divided into three groups (the black circles, red triangles and blue crosses) based 
on their proximity to the high-temperature limit, where k is expected to be the only figure of merit. 
The circled point with the arrow pointing to it corresponds to a simulation of the parameters for 
the last source in Table [ij which has k = 1.1 x 1CT 2 s" 1 (see main text). 
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The first two panels of Fig. [5] both display the equilibrium condensate number as a 
function of the phase-space flux k. Fig. |5](a) uses a log-log scale to highlight the behaviour 
for small and large values of k, while Fig. |5](b) uses a linear-linear scale to demonstrate that 
the black circles lying on a single line in Fig. [5|a) are not simply an artefact of plotting 
the results using a logarithmic scale. Finally, Fig. |5](c) displays the equilibrium condensate 
fraction function of k. 

Figures |5|a) and (b) demonstrate that it would be possible to produce atom lasers with 
respectable condensate numbers N > 10 5 (corresponding to atom laser fluxes of > 3 x 10 4 
atoms/s for the outcoupling rate 7 = 0.3 s" 1 chosen) by using replenishment sources that have 
a phase-space flux k > 10 -3 s . To determine if this is experimentally feasible, the properties 
of a range of experimental atomic sources are detailed in Table [T] and the corresponding values 
of the phase-space flux k calculated. 

The first source listed in Table [T] is the experiment of Ref. [6] that merged independently 
produced BECs in optical dipole traps that was discussed earlier. While this experiment is 
certainly not in the high-temperature limit, it has been included for comparison purposes. Of 
the remaining sources listed in Table [TJ most are many orders of magnitude away from being 
useful potential sources for a pumped atom laser (cf. Fig. [5]). The fluxes obtainable from 
these sources are insufficient to compensate for their higher temperatures as an increase of 
three-orders of magnitude in flux is necessary to compensate for an increase of a single order 
of magnitude in temperature (see ([9])). Only the last atomic source satisfies the requirement 
k > 10~ 3 s _1 . This experiment by Miiller et al. [7] is one of the sources in a dual atom 
interferometer designed for the precision measurement of accelerations and rotations |31j . A 
direct simulation has been performed for the parameters of this source, and the results are 
marked by a circle with an arrow pointing to it in Fig. [5j 

The steady-state condensate number for the source of Miiller et al. [7] is N = 5 x 10 5 
atoms, which would be a sufficiently large condensate to serve as a stable phase-reference 
for an atom laser, were it a pure BEC. However the equilibrium condensate fraction for this 
source is only 10% (see[5^c)). Previous theoretical work investigating the transfer of statistics 
from a trapped (quasi-)condensate to an atom laser found that using high-momentum kick 
Raman outcoupling such as that proposed in the scheme presented here can filter some 
of these fluctuations causing the atom laser to have a larger coherence length than the 
condensate from which it was produced [52] . 
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Atomic flux 


Temperature 


Phase-space flux 




Atomic source 


$ 


T 


K 




BECs in dipole traps 


10 5 s" 1 (55 mHz) a& 


< 1/iK 


> 1.9 x 10" 3 s" 1 


l — i 


2D+-MOT 


9 x 10 9 s" 1 


38 mK cd 


3 x 10~ 12 s" 1 


125 


2D+-MOT 


2 x 10 10 s" 1 


42 mK cd 


5 x 10~ 12 s" 1 


126 


MM-MOT 


10 9 s" 1 


61/uK c 


8 x 10" 5 s" 1 


[27 


LVIS 


5 x 10 9 s" 1 


25 mK cd 


6 x 10~ 12 s" 1 


[28] 


Zeeman slower 


3.2 x 10 12 s" 1 


32 mK cd 


2 x 10~ 9 s" 1 


129 

l S 


Magnetic guide loaded 


7 x 10 9 s" 1 


400//K c 


2 x 10~ 6 s 1 


8 


from 3D-MOT 










3D-MOT loaded from 


2 x 10 10 s" 1 (0.5 Hz) a 


500/xK 


3 x 10~ 6 s" 1 


m 


Zeeman slower 










3D-MOT loaded from 


3 x 10 s s" 1 (3 Hz) a 


8/xK 


1.1 x 10" 2 s" 1 




2D+-MOT 











a This source is pulsed, and the flux is the mean flux over one cycle with the repetition rate listed in 
parentheses. 

ft This repetition rate it too low for this source to be useful (see main text). It is listed for purposes of 
comparison only. 

c In keeping with the best-case scenario investigation being performed, this temperature assumes that the 

mean velocity of the atoms can be reduced to zero without affecting the distribution. This could be achieved, 

for example, by firing the source vertically below the main pumped atom laser experiment and taking the 

atoms from the mean turning point. 
d The dominant contribution to this temperature is the spread in the longitudinal velocities of the atoms. 

TABLE I: Relevant properties of selected experimental cold atomic sources. The phase-space flux 
k is evaluated from the atomic flux and temperature values listed using ([9]). 

It is not possible to investigate the transfer of statistics from the trapped component 
to the atom laser within the present model due to the simplifying assumption that it is 
only the condensate mode that is outcoupled to form the atom laser. A more detailed three- 
dimensional model taking into account the full spatial dependence of the Raman outcoupling 
process would be necessary to fully determine the feasibility of using an atomic source such 
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as that described by Miiller et al. [7J in the production of a truly continuous pumped atom 
laser. 

IV. CONCLUSIONS AND OUTLOOK 

We have investigated the feasibility of producing a continuously pumped atom laser 
fed by evaporation and replenishment with a thermal source. The method has been to 
investigate the best-case scenario in which the replenishment process introduces no heating 
to the trapped thermal component beyond that due to bringing the replenishing atoms into 
contact with the thermal cloud. With these caveats in mind, the results are promising: using 
an existing experimental source [7j it appears possible to produce steady-state condensates 
with large atom number (~ 5 x 10 5 atoms) using the scheme presented in Fig. [I] Could the 
atomic flux of this source be increased by an order of magnitude, the condensate number 
produced by this scheme could be pushed to 5 x 10 6 atoms. 

Ultimately it is not the size of the condensate that we are interested in, but the resulting 
flux of the atom laser and its coherence length. The former goal will certainly be improved 
by larger equilibrium condensate sizes, however the atom laser itself will be useless unless 
its coherence length is sufficiently large compared to MOT thermal sources to offset the 
reduced atom flux. It is by no means clear from the present work whether an atom laser 
produced using the atomic source of Miiller et al. [7] will have a coherence length that is 
larger than the size of the condensate from which it was outcoupled. To investigate this 
point it will be necessary to not make the assumption of ergodicity and to include the full 
spatial dependence of the Raman outcoupler. This could be achieved by making use of a 
finite temperature classical field theory [31] . 

One thing is certainly clear from the results presented here: the replenishment source for 
a collisionally-pumped atom laser must be close to degeneracy. It is simply not realistic to 
compensate for higher temperatures with a sufficiently increased flux. 

Acknowledgments 

This research was supported by the Australian Research Council Centre of Excellence for 
Quantum- Atom Optics (ACQ AO). We acknowledge the use of CPU time at the National 



19 



Computational Infrastructure National Facility. 



Appendix: Details of the Quantum Kinetic Theory model 

One of the difficulties involved in solving the kinetic model is that the energy range that 
the problem is defined over changes in time. The maximum energy is simply the energy 
of the evaporative cut-off £ cut , while the minimum energy is the chemical potential of the 
condensate p(t). A discretisation of the energy dimension over the range [0,e cut ] will suffer 
from problems accurately representing the lower-end of the distribution where the minimum 
energy of the thermal atoms is varying. An alternative is write the problem in terms of a 
shifted energy coordinate e = e — p(t) so that the minimum energy of the system is now fixed 
[12] ■ Of course the same problem now exists at the upper end of the energy range where 
the maximum energy e max = e cut — pit) is now time-dependent. However, at equilibrium 
there will be significantly fewer thermal atoms at the evaporation cut-off than there will be 
near the condensate (this is illustrated in Fig. [3^ a)). This choice will then result in smaller 
numerical errors than the alternative. 

Written in terms of the shifted energy variable e, the contribution due to the replenish- 
ment is 



where a bar over a function is used to indicate that it is defined in terms of the shifted 
energy coordinate. The original form of this term is given by (J3]). 

1. Density of states 

Although the evolution equations for the kinetic model ([I]) are written in terms of the 
product of the density of states p(e, t) and the energy distribution function g(e, t), it will be 
necessary to separately determine the energy distribution function to evaluate the collisional 
contributions given in the following section. To extract the energy distribution function it is 
necessary to have an explicit expression for the density of states for the thermal cloud. This 
density of states is not simply the same as that for a harmonic trap as the thermal modes 
will experience a mean-field repulsion due to the condensate mode. The effective potential 



d(p(e,t)g(e,t)) 
dt 




(A.I) 



replenishment 
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experienced by the thermal atoms is 



V eS (r, t) = Krap(r-) + 2gn c {r, t), (A.2) 

where VJ; rap (r) is the potential due to the magnetic trap, g = Air h 2 a/m, a is the s-wave 
scattering length and n c (r,t) is the condensate density which was assumed to follow a 
Thomas-Fermi distribution. Note that the '2' in the above expression is the full Hartree-Fock 
mean field experienced by the thermal atoms (see [T9J, Chapter 8] for further details) which 
is twice the mean-field repulsion experienced by condensate atoms. This is essentially due to 
the thermal atoms being distinguishable from the condensate atoms, while the condensate 
atoms are indistinguishable from one another. 



The density of states in the presence of the effective potential (A.2) is given by 



p(e, t) = J |^ 6 (e - V eS (r, t) - p 2 /2m) . (A.3) 

The integrals are performed in [TJ] giving the following result in terms of the shifted energy 
coordinate (Eqs. 49 and 50 in [12]) 



m*) = ^=[I-(e) + I + (e)], (A.4) 



where the functions I±(e) are 



3 o x=\/2u/h 

,_. u_x cl-U-X a_ v 

I -(e) = In [x + m_ 

v ' 4 8 8 V ; 



_ v?,x a + u + x a\ fx 
I+(e) = 1 h — arcsm 1 



v/ max{0,— a— } 
x= s Jaj r 



(A.5) 
(A.6) 

with a± = 2(e ± fi)/tujj, and u± = \/a± =F x 2 . Note that there is a minor typo in Bijlsma 
et al. [12], the lower limit of 1 -(e) is given as x = A/max{0, a_}, while it should read 
A/max{0, —a-} as in ( A.5[ ). 



x 



2. Collision and energy-redistribution in Quantum Kinetic Theory 

A full derivation of the forms of the collision and energy-redistribution terms of the kinetic 
model is given in [121 12T] . 
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The contribution due to thermal-thermal collisions is given in Eq. 26 of [12] and has the 
form 



d(p(si,t)g(e 1 ,t)) 



at 



3 2 

m g 
2vrW 



thermal-thermal 



de 2 / de 3 / de A p(e min ,t) 



(A.7) 



x S(ei + e 2 - £3 - £4) 

x [(1 + + 2 )»3»4 ~ 9x920- + 93)^ + 94)}, 
where £ m i n is the minimum of the £j, and g i = g(Ei,t). 

The contribution due to thermal-condensate collisions is given by Eq. 53 and Eq. 58-60 
of H21 and has the form 



9(p(si,t)g(ei,t)) 



dt 



thermal-condensate 



d&2 / des / cfe 4 5(e 2 — £3 — £4) 



3 2 
m g 

2ir 3 h 7 



x [5(e 1 - e 2 ) - S(ei ~ e 3 ) - 5(e 1 - e 4 )] 

X [(1 + 92)9394 ~ #2(1 + + 94)} 

dr n c (r, t), 

U eS (r,t)<U- 



(A.8) 



where U. 



[£3 + £4 



£3 — £3^4 + £4 



and U e {j{r,t) = U e g(r,t) — p{t). The corre- 



sponding contribution to the evolution of the condensate number is simply 



dN 



dt 



thermal-condensate 



- / de 



d(p(e,t)g(e,t)) 



dt 



(A.9) 



thermal-condensate 

Finally, the contribution due to energy redistribution is (Eqs. 32 and 52 in [12]) 



9(p(si,t)g(£i,t)) 



dt 

where p w is the weighted density of states 

2 



redistribution 



d(p w g) 
de 



(A.10) 



(A.ll) 



where the functions I±{e) are given in (A. 5) and (A. 6). 



3. Three-body loss in Quantum Kinetic Theory 

The dominant density-dependent loss process in Bose-Einstein condensates is three-body 
loss [35l [36] . Three-body loss (or three-body recombination) is the process in which three 
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atoms collide forming a bound dimer with the third necessary to ensure both energy and 
momentum conservation. The binding energy is sufficient to give the products of a three- 
body recombination process sufficient kinetic energy to rapidly escape the trap. Three-body 
loss is then well-described by the master equation term 

f 



dp 
dt 



L, / dxV 



(A.12) 



3-body loss 

where V[c]p = cp& — \{&cp + p&c) is the usual decoherence superoperator, and L 3 = 
5.8 x lCr 30 cm 6 s -1 [35] is the three-body recombination loss rate constant. This equation, 
first derived by [37] has the familiar form of a decoherence superoperator with the state 
undergoing loss as the argument. 

The loss rate of atoms from the system due to three-body loss is readily obtained from 



(A.12) as 



dN 
~dt 



3-body loss 



Tr { / dr&(r)V(r] 



dp 
dt 





f~* J 


3-body loss 





L 3 J dr{&(r) 3 V(r) 3 ). (A.13) 



To separate the contributions to (A.13) due to the thermal and condensed components, 
we use a broken symmetry approach. We write the annihilation operator ^ in terms of its 
mean value \& = (Jf?} and the fluctuation operator 5^/ = ^ — ^ and substitute this into 



(A.13). The fluctuation operator defined here includes thermal fluctuations, which cannot 
be considered to be small. Higher powers of 5^ can therefore not be neglected. However, 
thermal fluctuations have no well-defined phase relationship to one another or to the con- 
densate. Expectation values containing an unequal number of creation and annihilation 
deviation operators such as can therefore be assumed to be zero. 



Performing the substitution described, (A.13) becomes 



dN 
~dt 



3-body loss 



L 3 / dr |K(r)] 3 +9[n c (r)] 2 (5&(r)5i?(r)) 



+ 



9n c (r)(<J* t (r) 2 (J*(r) 2 > + (^(r) 3 ^^) 3 )}, 



(A.14) 



where n c (r) = \^(r)\ 2 is the condensate density. 

The non-condensate density is given by n^(r) = (SW (r)5^(r)V As thermal states are 
Gaussian, the higher-order expectation values in the previous expression may be simplified 
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by the application of Wick's theorem [38] giving 



(5¥(r) 2 5$(r) 2 ) = 2[n T (r)] 2 , 
(6&(r) 3 6*{rf) = 6[n T (r)f. 



(A.15) 
(A.16) 



Substituting these expressions back into (A. 14) yields 



^ =-Ls f dr K(r)] 3 + 9[n c (r)} 2 n T (r) + 18n c (r)[n T (r)] 2 + 6[n T (r)] 3 . (A.17) 

dt 3-body loss J 

The evaluation of this loss rate requires the evaluation of the condensate and thermal 
densities. The condensate density n c (r) is fully determined by the condensate occupation 
No(t) within the Thomas- Fermi approximation that has already been made elsewhere in 



the derivation of the kinetic model. The first term of (A.17) only involves the condensate 
density and may be evaluated analytically 

dN T 15 4 / 5 / 



, 12/5 

\ 9/5 



dt 



168vr 2 \hy/a 



(A.18) 



The remaining terms of (A.17) require an expression for the thermal density n^(r), which 



can be obtained from the energy distribution function g(e) and the density of states p(s). 
The total number of thermal atoms Nt can be written as 



N 7 



dep(e)g(e), 



(A.19) 



where the density of states is defined by (A. 3). Substituting this into (A.19) and rearranging 
the order of integrals gives 



N'T 



dr / dep(e,r)g(e), 



(A.20) 



where we have defined 



/™3/2 
dp 5 (e - V ee (r, t) - p 2 /2m) = -^-^ y/e-V^r). (A.21) 



The thermal density can be identified from (A.20) 



n T {r) = J dep(£,r)g(e). 



(A.22) 



The remaining terms of (A.17) can now be expressed in terms of the energy distribution 



function g{e) and the density of states p(e) by substituting (A.22) for one of the factors of 
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riT{r) in each term 

-L 3 J dr9[n c (r)] 2 n T (r) = -L 3 J de g(e) J dr9p(e,r)[n c (r)] 2 , (A.23) 

-L 3 J drl8n c {r)[n T {r)] 2 = -L 3 J de g(e) J dr 18 p(e,r)n c (r)n T (r), (A.24) 

-L 3 J drQ[n T (r)f = -L 3 J de g(e) J dr Qp(e,r)[n T (r)} 2 . (A.25) 

From these expressions the rate of loss of atoms of energy e from the distribution can be 
identified 

= -L 3 J drp(e,r)g(e){3[n c (r)] 2 + 12n c (r)n T (r) + Q[n T (r)] 2 }, 

3-body loss 

(A.26) 



d( P (e)g(e)) 



dt 



where the contributions due to the terms involving only one or two thermal atoms have 
been multiplied by 1/3 and 2/3 respectively to share appropriately the total loss. The 
corresponding term for the condensate number evolution is 

15 4 / 5 / muj \ 12/5 



dN 



9/5 



dt 3-body loss L3 1687T 2 \hy/a) N ° (A 2?) 

-L 3 J de J drp(e,r)g(e) {6[n c (r)] 2 + Qn c (r)n T (r)} , 

where the contributions due to the terms involving only one or two condensate atoms have 
been multiplied by 1/3 and 2/3 respectively. 
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